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Abstract 

In a recent publication Zheng, Vary, and Barrett reproduced the negative 
quadrupole moment of 6 Li and the low-lying positive-parity states of 5 He 
by using a no-core shell model. In this Comment we question the meaning 
of these results by pointing out that the model used is inadequate for the 
reproduction of these properties. 
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Recently, Zheng, Vary and Barrett have published a number of no-core shell-model calcu- 
lations for light nuclei in a series of papers; see Q and references therein. In their most recent 
article |]J they presented results for two long-standing problems of light nuclei. Namely, they 
reproduced the negative quadrupole moment of 6 Li and found the previously predicted low- 
lying positive-parity states of 5 He. In this Comment we argue that these results cannot be 
considered the solutions to the problems in question since the model is not realistic enough 
to cope with these problems. We first discuss the quadrupole moment of 6 Li and then the 
positive-parity states of 5 He. 

The experimental value of the 6 Li quadrupole moment is —0.083 e-fm 2 0. There have 
been numerous theoretical calculations attempting to reproduce this negative value, without 
success. They include variational three-body calculations |§, the hyperspherical harmonics 
expansion method [f|], the Faddeev approach [0], and a six-body a + p + n three-cluster 
model |J. All these calculations gave positive values for Q, in the 0.2-0.6 e-fm 2 range. In 
Ref. [0 it was claimed that a three-cluster model can reproduce the negative quadrupole 
moment of 6 Li if the full six-body antisymmetrization is properly taken into account. This 
claim was, however, disputed ||, and the result was shown to arise from a restriction on the 
model space (see also ||). These macroscopic and microscopic a+p+n calculations revealed 
that the value of the 6 Li quadrupole moment results from a delicate balance between the 
contributions coming from configurations of different angular momenta in the p-n and at- 
(pn) relative motions @||. This indicates that, for Q to be reliable, the p-n and a-(pn) 
dynamics must be described correctly. 

In Ref. [0J the model produces Q = —0.116 e-fm 2 . It would be the first time that a 
realistic microscopic model reproduced the correct negative sign of Q. However, in Table I 
of Ref. m we can see that the ground-state binding energy of 6 Li is incompatible with those 
of the deuteron and of the a particle. As a consequence, the model 6 Li seems to be unstable 
against the a + d breakup by 0.21 MeV, in sharp contrast with reality. (The ground state 
of 6 Li is below the a + d threshold by 1.475 MeV 0.) In fact, the separation energy should 
be calculated from binding energies obtained in fully consistent state spaces, which is not 
exactly so in Ref. 0. The state space used for 6 Li is restricted by 2hu more than those of a 
and d, thus the breakup energy implied by the 6 Li wave function may differ from 0.21 MeV. 
Nevertheless, it looks likely that such a major discrepancy cannot result from a such a minor 
mismatch between the state spaces. To reproduce the sort of loose intercluster binding that 
is characteristic of the a-d relative motion is a major challenge to any shell model ||. 

Since the 6 Li wave function contains several configurations, each belonging to a spe- 
cific separation energy, it would be difficult to demonstrate the strong dependence of the 
quadrupole moment on all these separation energies. Just to give an indication for the mag- 
nitude of such a dependence, in Fig. 1 we show the quadrupole moment of the ground state 
of 7 Li in ana + t cluster model, using the Minnesota effective nucleon-nucleon force [|K]] . To 
mimic the situation in Ref. the exchange mixture parameter of the interaction is varied 
such that the energy of 7 Li changes, while the binding energies of a and t remain fixed. 
We show Q as a function of the a+t separation energy. One can see that the quadrupole 
moment is very sensitive to this quantity. The 6 Li quadrupole moment is likely to behave 
in a similar way, so that an increase of the a + d separation energy by some 1.68 MeV could 
considerably shift it toward positive values. 

In summary, a theoretical value for the 6 Li quadrupole moment can only be accepted as a 
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physically meaningful value if the model produces the correct separation energies, especially 
for the most important a + d partition. Without this, a reproduction of the experimental 
value can only be fortuitous. 

The problem with the low-lying positive-parity states of 5 He in Ref. is that the model 
cannot handle the asymptotics of the wave function correctly. In Ref. |1[ the zero point of 
the energy scale is at the total A-body disintegration threshold of the A-nucleon system. So, 
any state which is below this A-body threshold has a negative energy, which would imply 



an exponentially damped bound-state v4-body asymptotic wave function [11 



^ A ~ exp(-k A p A ) for p A -> oo. (1) 

Here p A is the hyperradius, p\ = Y^irf, where rj are the one-particle position vectors, and 
k A = (2rri]yE A /h ) 1 ' 2 , where E A is the binding energy of the A-body system relative to the 
A-body disintegration threshold, and is the nucleon mass. Such a ^ A can really be 
expanded in terms of square-integrable functions, like in a shell model. However, the Eq. (1) 
boundary condition is correct, describing a physical situation, only if there is no break-up 
channel below E A . If there is a two-body (A = B + C) break-up channel below E A , then 
the correct boundary condition is [12| 



^> A ~ exp(— k A p A ) + [xexp(— ikr) + y exp(i/cr)] 

for p A , r — > oo. (2) 

Here r is the distance between the fragments B and C, the functions $ are the internal 
states of B and C with binding energies Eb and Ec, and k = [2p(E A — Eb — Ec) /h 2 } 1 ^ 2 . A 
scattering "state", of energy E A , that obeys Eq. (2) will be regarded as a (resonant) state 
of the nucleus if S = —y/x (the "S-matrix") has a pole at the complex energy E A — iT/2, 
where T is the total width. 

Square-integrable bases, like that of Ref. 0, are obviously unable to observe the bound- 
ary condition in Eq. (2), so whatever they predict for states above break-up thresholds is to 
be taken with reservation. This applies to the low-lying positive-parity states of 5 He, which 
are above the a + n threshold. 

We have performed a search for such states in a cluster model whose wave function is 

* = Y.[M[{^ n )sxT{p a n)} JM } 

S,L 

+ A{[(<S> d &)sxf(p dt )} JM }}, (3) 

where A is the intercluster antisymmetrizer, the cluster internal states $ are translation 
invariant Os harmonic oscillator shell-model states for the a particle, deuteron, and triton, 
the vectors p are the intercluster Jacobi coordinates, and [...] denotes angular momentum 
coupling. In the sum over S and L all possible configurations are included. The intercluster 
relative-motion functions \ have the correct asymptotics. To avoid any ambiguity in the 
recognition of a resonance in the phase shift 5 = \ arg S, we searched for complex-energy 



poles of the S-matrix directly. Both an analytic continuation method [13| and the complex 



scaling method |14| were used. The 3/2 and 1/2 states were found, but the next level was 



the 3/2 + state at 16 MeV excitation energy. No sign for any low-lying l/2 + , 3/2 + , or 5/2" 
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states was found. The inclusion of a few monopole breathing excitations of a did not change 
the situation either. This rules out even the exotic possibility that the low-lying positive- 



parity states are Pauli resonances [15], since with the departure from the single oscillator 
description of the a cluster, the configurations that might produce Pauli resonances get 
automatically included. 

To produce a low-energy l/2 + state artificially, we made the intercluster binding stronger, 
while keeping the cluster binding energies fixed, by changing a mixing parameter of the N-N 
interaction. Then we let the mixing parameter tend toward its physical value, and followed 
the position of the l/2 + pole. The pole moved rapidly towards higher energies and, for 
example, it was found at 19 MeV excitation energy with a width of 34 MeV while the 
mixing parameter was still highly non-physical; indeed, the same parameter value produced 
deeply bound 3/2~ and 1/2 - states. Further change of the mixing parameter in the direction 
of its correct value pushed this l/2 + state to even higher energies with larger widths. 

All in all, in a model that handles the asymptotics correctly, the low-lying positive-parity 
states of 5 He do not show up. Although the basis used by Zheng et al. |l] is probably more 
flexible than ours to describe the correlated short-range motion of the nucleons, it is very 
difficult to imagine that an improvement of our model in this direction would bring down 
high-lying positive-parity states from the upper region of the continuum. Nevertheless, it 
is still to be proven whether or not the non-physical boundary condition results in the 
appearance of these states. The best way to check this would be to supplement the wave 
function of Ref. [|J by an a + n cluster term which describes the correct asymptotics. This 
could be done, for example, in a cluster-configuration shell- model jTjJ. An indirect indication 
as to whether these states are real or spurious could be obtained more simply by examining 
the stability of their energies against changes in the size of the square-integrable basis. 
Resonant states produced by a square-integrable basis should show stability |T7| ]. 

Just to give an example that the incorrect boundary condition can incur spurious states, 
here we show the case of 8 Li. This nucleus is described in a three-cluster a + t + n model 
with the basis containing a number of different angular momenta ||18|| . Expanding all inter- 
cluster relative-motion functions in terms of square-integrable functions, one gets, from the 
diagonalization of the Hamiltonian, two 2 + states below the a + t + n threshold, at binding 
energies 4.2 MeV and 1.1 MeV (relative to the three-cluster threshold), respectively. Experi- 
ment has only produced one 2 + state in this region, the ground state. The square-integrable 
basis is seemingly adequate, because we are below the three-cluster threshold. However, the 
1.1 MeV-energy state is above the 7 Li+n two-body break-up threshold, which means that 
the correct boundary condition has to contain a 7 Li+n scattering term. Supplementing our 
wave function by such a term, the 1.1 MeV 2 + state disappears immediately, showing that 
this state was an artifact brought about by the incorrect boundary condition. 

By this analogy, we suggest that the low-lying positive-parity states of 5 He could also 
disappear if the proper boundary condition were taken into account. Be that as it may, we 
cannot claim that this would disprove the existence of the states in question. But certainly, 
before further efforts are spent on understanding the nature of these states, the empirical 
evidence for their existence should be reconsidered. 

In conclusion, the reproduction of the negative quadrupole moment of 6 Li and the low- 
lying positive-parity states of 5 He by Zheng et al. in a shell model cannot be regarded as 
well-founded because the aspects of few-body dynamics underlying these effects are treated 
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improperly in that model. 
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FIGURES 

FIG. 1. The quadrupole moment of the ground state of 7 Li as a function of its a + t separation 
energy in an a + t cluster model. The experimental separation energy is 2.47 MeV, and the 
quadrupole moment is —4.06 e-fm 2 [p|. 
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